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Wigner's quasi-probability distribution function in phase-space is a special (Weyl) rep- 
resentation of the density matrix. It has been useful in describing quantum transport 
in quantum optics; nuclear physics; decoherence (eg, quantum computing); quantum 
chaos; "Welcher Weg" discussions; semiclassical limits. It is also of importance in signal 
processing. Nevertheless, a remarkable aspect of its internal logic, pioneered by Moyal, 
has only emerged in the last quarter-century: It furnishes a third, alternative, formu- 
lation of Quantum Mechanics, independent of the conventional Hilbert Space, or Path 
Integral formulations. In this logically complete and self-standing formulation, one need 
not choose sides — coordinate or momentum space. It works in full phase-space, accom- 
modating the uncertainty principle. This is an introductory overview of the formulation 
with simple illustrations. 

1. Introduction 

There are three logically autonomous alternative paths to quantization. The first 
is the standard one utilizing operators in Hilbert space, developed by Heisenberg, 
Schrodinger, Dirac, and others in the Twenties of the past century. The second one 
relies on Path Integrals, and was conceived by Dirac^ and constructed by Feynman. 

The last one (the bronze medal!) is the Phase-Space formulation, based on 
Wigner's (1932) quasi-distribution function^ and Weyl's (1927) correspondence'^ 
between quantum-mechanical operators and ordinary c-number phase-space func- 
tions. The crucial composition structure of these functions, which relies on the 
★-product, was fully understood by Groenewold (1946)^, who, together with Moyal 
(1949)"", pulled the entire formulation together. Still, insights on interpretation and 
a full appreciation of its conceptual autonomy took some time to mature with the 
work of Takabayasi^, Baker^, Fairlie ®, and others. (A brief guide to some landmark 
papers is provided in the Appendix.) 

This complete formulation is based on the Wigner Function (WF), which is a 
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quasi-probability distribution function in phase-space, 

f{x,p) ^^Jdy r{x - ^y) e-'yp^ix + ^y). (1) 

It is a special representation of the density matrix (in the Weyl correspondence, 
as detailed in Section 8). Alternatively, it is a generating function for all spatial 
autocorrelation functions of a given quantum- mechanical wavefunction ip{x). 

There are several outstanding reviews on the subject'^'"'^'^'"'^"'^. Nevertheless, the 
central conceit of this review is that the above input wavefunctions may ultimately 
be forfeited, since the WFs are determined, in principle, as the solutions of suitable 
(celebrated) functional equations. Connections to the Hilbert space operator formu- 
lation of quantum mechanics may thus be ignored, in principle — even though they 
are provided in Section 8 for pedagogy and confirmation of the formulations' equiv- 
alence. One might then envision an imaginary planet on which this formulation 
of Quantum Mechanics preceded the conventional formulation, and its own tech- 
niques and methods arose independently, perhaps out of generalizations of classical 
mechanics and statistical mechanics. 

It is not only wavefunctions that are missing in this formulation. Beyond an 
all-important (noncommutative. associative, pseudo differential) operation, the 
product, which encodes the entire quantum mechanical action, there are no opera- 
tors. Observables and transition amplitudes are phase-space integrals of c-number 
functions (which compose through the ★-product), weighted by the WF, as in sta- 
tistical mechanics. Consequently, even though the WF is not positive-semidefinite 
(it can, and usually does go negative in parts of phase-space), the computation of 
observables and the associated concepts are evocative of classical probability theory. 

This formulation of Quantum Mechanics is useful in describing quantum trans- 
port processes in phase space, of importance in quantum optics, nuclear physics, 
condensed matter, and the study of semiclassical limits of mesoscopic systems^^ 
and the transition to classical statistic;al mechanics^^. It is the natural language 
to study quantum chaos and decoherence^^ (of utility in, eg, quantum computing), 
and provides intuition in QM interference such as in "Welcher Weg" problems^^, 
probability flows as negative probability backflows^^, and measurements of atomic 
systems""^^. The intriguing mathematical structure of the formulation is of relevance 
to Lie Algebras^*, and has recently been retrofitted into M-theory advances linked 
to noncommutative geometry and matrix models^", which apply spacetime un- 
certainty principles'^ reliant on the ★-product. (Transverse spatial dimensions act 
formally as momenta, and, analogously to quantum mechanics, their uncertainty is 
increased or decreased inversely to the uncertainty of a given direction.) 

As a significant aside, the WF has extensive practical applications in signal 
processing and engineering (time- frequency analysis), since time and energy (fre- 
quency) constitute a pair of Fourier-conjugate variables just like the y and p of phase 
space? 

"Thus, time Vcirying signals are best represented in a WF as time varying spectrograms, analogously 
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For simplicity, the formulation will be mostly illustrated for one coordinate 
and its conjugate momentum, but generalization to arbitrary-sized phase spaces is 
straightforward^^, including infinite-dimensional ones, namely scalar field theory^^'^^'^^'' 

2. The Wigner Function 

As already indicated, the quasi-probability measure in phase space is the WF, 

fix,p) = ^Jdy r{x - \y) e-'yp^{x + ^y). (2) 

It is obviously normalized, J dpdxf{x,p) = 1. In the classical limit, fi, — > 0, it would 
reduce to the probability density in coordinate space x, usually highly localized, 
multiplied by (5-functions in momentum: the classical limit is "spiky"! This ex- 
pression has more x — p symmetry than is apparent, as Fourier transformation to 
momentum-space wavefunctions yields a completely symmetric expression with the 
roles of x and p reversed, and, upon rescaling of the arguments x and p, a symmetric 
classical limit. 

It is also manifestly Real? It is also constrained by the Schwarz Inequality to be 
bounded, — f < f{x,p) < f • Again, this bound disappears in the spiky classical 
limit. 

p- or .T-projcction leads to marginal (bona-fide) probability densities: a spacclike 
shadow J dpf{x,p) = p{x), or else a momentum-space shadow / dxf{x,p) — <y{p), 
respectively. Either is a (bona-fide) probability density, being positive semidefinite. 
But neither can be conditioned on the other, as the uncertainty principle is fighting 
back: The WF f{x,p) itself can, and most often does go negative in some areas 
of phase-space^'^, as is illustrated below, a hallmark of QM interference in this 
language. (In fact, the only pure state WF which is non-negative is the Gaussian^^.) 

The counter-intuitive "negative probability" aspects of this quasi-probability 
distribution have been explored and interpretcd^'^''^*''^^, and negative probability 
flows amount to legitimate probability backflows in interesting settings^^. Never- 
theless, the WF for atomic systems can still be measured in the laboratory, albeit 
indirectly 

Smoothing / by a filter of size larger than H (eg, convolving with a phase-space 
Gaussian) results in a positive-semidefinite function, ie it may be thought to have 

been coarst-iu'd to a classical distribution^"^? 

to a music score, ic the changing distribution of frequencies is monitored in timc^^; even though 
the description is constrained and redundant, it gives an intuitive picture of the signal that a mere 
time profile or frequency spectrogram fails to convey. Applications abound in acoustics, speech 
analysis, seismic data analysis, internal combustion engine-knocking analysis, failing helicopter 
component vibrations, etc. 

''In one space dimension, by virtue of non-degeneracy, it turns out to be p-even, but this is not a 
property used here. 

'^This one is called the Husimi distribution^^, and sometimes information scientists examine it on 
account of its non-negative feature. Nevertheless, it comes with a heavy price, as it needs to be 
"dressed" back to the WF for all practical purposes when expectation values axe computed with 
it, ie it does not serve as an immediate quasi-probability distribution with no other measure^^. 
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Among real functions, the WFs comprise a rather small, highly constrained, 
set. When is a real function f{x,p) a bona-fide Wigner function of the form (2)? 
Evidently, when its Fourier transform "left-right" factorizes, 

fix,y) = j dp e'pyf{x,p) = gl{x - hy/2) gnix + hy/2), (3) 



le, 



-0, (4) 



d{x - hy/2) d{x + hy/2) 

so gL = gn from reality. 

Nevertheless, as indicated, the WF is a distribution function: it provides the 
integration measure in phase space to yield expectation values from phase-space 
c-number functions. Such functions are often classical quantities, but, in general, 
are associated to suitably ordered operators through Weyl's correspondence rule^. 
Given an operator ordered in this prescription, 

©(j,P) = -^2:^ / ^'^d.adxdp g{x,p)exp{iT{p-p) + ia{f- x)), (5) 

the corresponding phase-space function g{x,p) (the "classical kernel of the opera- 
tor" ) is obtained by 

p^p, y a;. (6) 
That operator's expectation value is then a "phase-space average" 

(©) = / dxdp f{x,p) g{x,p). (7) 



The classical kernel g{x,p) is often the unmodified classical expression, such as a 

2 

conventional hamiltonian, H = ^ + ic the transition from classical me- 

chanics is straightforward ; but it contains h when there are quantum-mechanical 
ordering ambiguities, such as for the classical kernel of the square of the angular 
momentum £ • £: this one contains 0{h'^) terms introduced by the Weyl ordering, 
beyond the classical expression L^. In such cases, the classical kernels may still be 
produced without direct consideration of operators. A more detailed discussion of 
the correspondence is provided in Section 8. 

In this sense, expectation values of the physical observables specified by classical 
kernels g{x,p) are computed through integration with the WF, in close analogy to 
classical probability theory, but for the non-positive-definiteness of the distribution 
function. This operation corresponds to tracing an operator with the density matrix 
(Sec 8). 

The negative feature of the WF is, in the last analysis, an asset, and not a liability, and provides 
an efficient description of "beats" 
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3. Solving for the Wigner Function 

Given a specification of observables, the next step is to find the relevant WF for a 
given hamiltonian. Can this be done without solving for the Schrodinger wavefunc- 
tions tp, i.e. not using Schrodinger's equation directly? The functional equations 
which / satisfies completely determine it. 

Firstly, its dynamical evolution is specified by Moyal's equation. This is the 
extension of Liouville's theorem of classical mechanics, for a classical hamiltonian 
H{x,p), namely dtf + {/, H} = 0, to quantum mechanics, in this language^'^: 

where the ★-product^ is 

-k = eridxdp-dpdx) _ (cf^ 

The right hand side of (8), dubbed the "Moyal Bracket" (the Weyl-correspondent 
of quantum commutators), is the essentially unique one-parameter (7i) associative 
deformation of the Poisson Brackets of classical mechanics^^ . Expansion in h around 
reveals that it consists of the Poisson Bracket corrected by terms 0{h). The equa- 
tion (8) also evokes Heisenberg's equation of motion for operators, except H and 
/ here are classical functions, and it is the ;*r-product which introduces noncommu- 
tativity. This language, thus, makes the link between quantum commutators and 
Poisson Brackets more transparent. 

Since the ★-product involves exponentials of derivative operators, it may be 
evaluated in practice through translation of function arguments ("Bopp shifts"), 

f{x,p)*g{x,p) = f ix + — dp, p - — dxj g{x,p). (10) 

The equivalent Fourier representation of the ★-product is^: 

f^9 = J dp'dp"dx'dx" f{x',p') gix",p") 

= exp(^^{p{x' -x")+p\x" -x)+p"{x-x'))y (11) 

An alternate integral representation of this product is^^ 



f-kg={h-n) ^ j dudvdwdz f{x+u,p+v) g{x+w,p+z) exp 



I {uz -vw)], (12) 



which readily displays noncommutativity and associativity. 

★-multiplication of c-mimber phase-space functions is in complete isomorphism 
to Hilbert-space operator multiplication^, 



(27r)2 / 



2l*B = „ / drdadxdp (a ★ b) exp(iT(p — p) + icF{l — x)). (13) 
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The cyclic phase-space trace is directly seen in the representation (12) to reduce to 
a plain product, if there is only one ★ involved, 



Moyal's equation is necessary, but does not suffice to specify the WF for a system. 
In the conventional formulation of quantum mechanics, systematic solution of time- 
dependent equations is usually predicated on the spectrum of stationary ones. Time- 
independent pure-state Wigner functions ★-commute with H, but clearly not every 
function ★-commuting with H can be a bona-fide WF (eg, any ★-function of H will 
★-commute with H). 

Static WFs obey more powerful functional ★-genvalue equations® (also see ^^), 



where E is the energy eigenvalue of = -^'V'- These amount to a complete 
characterization of the WFs'^^. 

Lemma: For real functions f{x,p), the Wigner form (2) for pure static eigenstates 
is equivalent to compliance with the -k-genvalue equations (15) (5R and 9 parts). 
Proof: 




(14) 




(15) 



H{x,p)-k f{x,p) = 




= E f{x,p). 



(16) 



Action of the effective differential operators on ij;* turns out to be null. 
Symmetrically, 
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= E f{x,p), (17) 

where the action on ip is now trivial. Conversely, the pair of ^-eigenvalue equations 
dictate, for f{x,p) = Jdy e-'^P/(x,y) , 



J dy e-'yp (-^{dy ±^ d.f + V{x±^y) - f{x,y) = 0. 



(18) 



Hence, Real solutions of (15) must be of the form 

f = Jdy e-'yp^*ix - ^y)i){x + |y)/27r, such that i^V^ = E^j. □ 

The eqs (15) lead to spectral properties for WFs^'"^®, as in the Hilbert space 
formulation. For instance, projective orthogonality of the it-genf unctions follows 
from associativity, which allows evaluation in two alternate groupings: 

fi.H*g = Ef f*g = Eg f*g. (19) 

Thus, for Eg ^ Ef, it is necessary that 

/*5 = 0. (20) 

Moreover, precluding degeneracy (which can be treated separately), choosing f = g 
above yields, 

f*H*f = Eff^f = H*f^f, (21) 
and hence / * / must be the stargenfunction in question, 

/*/<x/. (22) 
Pure state fs then ★-project onto their space. In general, it can be shown^'^^ that 

fa* fb=^ Sa.fi fa. (23) 

The normalization matters^: despite linearity of the equations, it prevents super- 
position of solutions. (Quantum mechanical interference works differently here, in 
comportance with density matrix formalism.) 

By virtue of (14), for different ★-genfunctions, the above dictates that 

dpdx fg = . (24) 

Consequently, unless there is zero overlap for all such WFs, at least one of the two 
must go negative someplace to offset the positive overlap^ — an illustration of the 
negative values' feature, quite far from being a liability. 
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Further note that integrating (15) yields the expectation of the energy, 

j H{x,p)fix,p) dxdp = E j f dxdp = E. (25) 

NB Likewise^ integrating the above projective condition yields 

f dxdp = ^ , (26) 

that is the overlap increases to a divergent result in the classical limit, as the WFs 
grow increasingly spiky. 

In classical (non-negative) probability distribution theory, expectation values 
of non-negative functions are likewise non-negative, and thus result in standard 
constraint inequalities for the constituent pieces of such functions, such as, e.g., 
moments of the variables. But it was just seen that for WFs which go negative for 
an arbitrary function g, need not > 0: this can be easily seen by choosing 

the support of g to lie mostly in those regions of phase-space where the WF / is 
negative. 

Still, such constraints are not lost for WFs. It turns out they are replaced by 

{9**9)>0, (27) 

and lead to the uncertainty principle^^. In Hilbert space operator formalism, this 
relation would correspond to the positivity of the norm. This expression is non- 
negative because it involves a real non-negative integrand for a pure state WF 
satisfying the above projective conditionf 

j dpdx{g*-kg)f = h j dxdp{g*-kg){f-kf) = h j dxdp{f-kg*)-k{g-kf) = h j dxdp\g-kf\'^ . 

(28) 

To produce Heisenberg's uncertainty relation, one only need choose 

g = a + bx + cp, (29) 

for arbitrary complex coefficients a, b, c. The resulting positive semi-definite quadratic 
form is then 



a*a+b*b{xi<x) + c*c{pikp) + {a*b+b*a){x) + {a*c+c*a){p)+c*b{p-kx) + b*c{xi^p) > 0, 
(30) 

''This discussion applies to proper WFs, eorresponding to pure states' density matrices. E.g., a 
sum of two WFs is not a pure state in general, and does not satisfy the condition (4). For such 
generalizations, the "impurity" is** J dxdp{f — hf^) > 0, where the inequality is only saturated 
into an equality for a pure state. For instance, for w s (/a + /i))/2 with /a * /i = 0, the impurity 
is nonvanishing, J dxdp{w — hw^) = 1/2. 

■^Similarly, if /i and /2 are pure state WFs, the transition probability between the respec- 
tive states is also non- negative^* , manifestly by the same argument^'' , namely, / dpdxfi fi = 
{2-KhfS dxdp |/i*/2p >0. 
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for any a, b, c. The eigenvalues of the corresponding matrix are then non-negative, 
and thus so must be its determinant. Given 

x-kx = x'^, P*P = P^, p-k X = px — ih/2, x-kp = px + ih/2, (31) 

and the usual 

{Axf ^ {{X - {xm, {Apf ^ Hp - (p))2), (32) 

this condition on the 3x3 matrix determinant amounts to 

{Axf {Apf > ny4+[{{x-{x))ip-{p)))y, (33) 

and hence 

A.X- Ap > h/2. (34) 
The h entered into the moments' constraint through the action of the ★-product. □ 



4. Illustration: the Harmonic Oscillator 

To illustrate the formalism on a simple prototype problem, one may look at the 
harmonic oscillator. In the spirit of this picture, one can, in fact, eschew solving 
the Schrodinger problem and plugging the wavefunctions into (2); instead, one may 
solve (15) directly for H = {p^ + x^)/2 (with m = 1, a; = 1), 



((^ + '^d,r + {p- '^d^f - 2E^ f{x,p) = 0. 



(35) 



For this hamiltonian, the equation has collapsed to two simple PDEs. The first one, 
the Imaginary part, 

{xdp-pd,)f = 0, (36) 

restricts / to depend on only one variable, the scalar in phase space, z = 4H/h = 
2{x'^ +p^)/h. Thus the second one, the Real part, is a simple ODE, 

t-zdl-d,-^\f{z) = Q. (37) 



,4 ' ^ 

Setting /(z) = exp(— z/2)L(z) yields Laguerre's equation 

E 1 



zdl + {l-z)d, + --^-] L{z) = 0. (38) 



It is solved by Laguerre polynomials. 



L„ = '^^^. (39) 
n! 
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for n = E/h — 1/2 = 0, 1, 2, .., so the ★-gen-Wigner-functions are 

Lo = 1, Li = l-4:H/h, L2=8H^/h^ -8H/h+l,... (40) 

But for the Gaussian ground state, they aU have zeros and go negative. These 
functions become spiky in the classical limit /i — s- 0; e.g., the ground state Gaussian 
/o goes to a (5-function. 

Their sum provides a resolution of the identity^, 

Y.fn = l/h. (41) 

n 

(For the rest of this section, for algebraic simplicity, set h — 1.) 

Oscillator Wigner Function, n=3 




Fig. 1. WF for the 3rd excited state. Note the negative values. 

Dirac's hamiltonian factorization method for the alternate algebraic solution of 
the same problem carries through intact, with *- multiplication supplanting operator 
multiplication. That is to say, 

H =^{x-ip)*ix + ip) + ^ . (42) 
10 
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This motivates definition of raising and lowering functions (not operators) 

a=-^{x + ip), a^ = -^{x-ip), (43) 

where 

a-ka^-a^*a = l . (44) 
The annihilation ones ★-annihilate the ★-Fock vacuum, 

a^/o = ^(a; + ip)*e-(^'+f') =0 . (45) 
v2 

Thus, the associativity of the ★-product permits the customary ladder spectrum 
generation^^. The ★-genstates for if ★ / = / ★ iJ are 

/„ cx (at^)" /o (★«)" . (46) 

They are manifestly real, like the Gaussian ground state, and left-right symmetric; 
it is easy to see they are ★-orthogonal for different eigenvalues. Likewise, they can 
be seen by the evident algebraic normal ordering to project to themselves, since the 
Gaussian ground state does, /o*/o = fo/h. The corresponding coherent state WFs 
^^'^^ are likewise very analogous to the conventional formulation. 

This type of analysis carries over well to a broader class of problems^^ with "es- 
sentially isospectral" pairs of partner potentials, connected with each other through 
Darboux transformations relying on Witten superpotentials W (cf the Poschl- Teller 
potential). It closely parallels the standard differential equation structure of the 
recursive technique. That is, the pairs of related potentials and corresponding ★- 
gcnstatc Wigncr functions arc constructed recursively'*^ through ladder operations 
analogous to the algebraic method outlined for the oscillator. 

Beyond such recursive potentials, examples of further simple systems where the 
★-genvalue equations can be solved on first principles include the linear potential 
41,36,42^ the exponential interaction Liouville potentials, and their supcrsymmctric 
Morse generalizations^^ (also see ^'^). Further systems may be handled through the 
Chebyshev-polynomial numerical techniques of ref 

First principles phase-space solution of the Hydrogen atom is less than straight- 
forward and complete. The reader is referred to ii>45,46 j^j. significant partial results. 

5. Time Evolution 

Moyal's equation (8) is formally solved by virtue of associative combinatoric opera- 
tions completely analogous to Hilbert space Quantum Mechanics, through definition 
of a ★-unitary evolution operator, a "★-exponential"^"'^ 

U4x,p;t) = e^/'^ ^ 1 + (it/h)H{x,p) + Wlff H + ^^^H H i. H + (47) 
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for arbitrary hamiltonians. The solution to Moyal's equation, given the WF at 
i = 0, then, is 

f{x,p;t) = U-\x, p;t)*f{x, p; 0) * U^{x, p; t) . (48) 

For the variables x and p, the evolution equations collapse to classical trajecto- 
ries, 

dx X -k H — II -k X ,^ 

-dt= -n = ^^^ = ^' ^''^ 

dp p-k H — H -k p 



dt ih 



= -d^H = -X , (50) 



where the concluding member of these two equations hold for the oscillator only. 
Thus, for the oscillator, 

x{t) = xcost + psint, (51) 



p{t) = pcost — xsint. (52) 

As a consequence, for the oscillator, the functional form of the Wigner function 
is preserved along classical phase-space trajectories^, 

f{x,p; t) = f {xcost — psint, pcost + xsint; 0). (53) 

Any oscillator WF configuration rotates uniformly on the phase plane around 

the origin, essentially classically, even though it provides a complete quantum me- 
chanical description4'4^'2'25.27,48 

Naturally, this rigid rotation in phase-space preserves areas, and thus automat- 
ically illustrates the uncertainty principle. By contrast, in general, in the conven- 
tional formulation of quantum mechanics, this result is deprived of intuitive import, 
or, at the very least, simplicity: upon integration in x (or p) to yield usual proba- 
bility densities, the rotation induces apparent complicated shape variations of the 
oscillating probability density profile, such as wavepacket spreading (as evident in 
the shadow projections on the x and p axes of Fig 2). 

Only when (as is the case for coherent states^^) a Wigner function configuration 
has an additional axial x — p symmetry around its own center, will it possess an 
invariant profile upon this rotation, and hence a shape-invariant oscillating proba- 
bility density^^. 

In Dirac's interaction representation, a more complicated interaction hamilto- 
nian superposed on the oscillator one, leads to shape changes of the WF config- 
urations placed on the above "turntable", and serves to generalize to scalar field 
theory^^. A more elaborate discussion of propagators in found in ^O'^^. 
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X 

Fig. 2. Time evolution of generic WF configurations driven by an oscillator hamiltonian. The 
t-arrow indicates the rotation sense of x and p, and so, for fixed x and p axes, the WF shoebox 
configurations rotate rigidly in the opposite direction, clockwise. (The sharp angles of the WFs 
in the illustration are actually unphysical, and were only chosen to monitor their "spreading 
wavepacket" projections more conspicuously.) These x and p-projections (shadows) are meant 
to be intensity profiles on those axes, but are expanded on the plane to aid visualization. The 
circular figure represents a coherent state, which projects on either axis identically at all times, 
thus without shape alteration of its wavepacket through time evolution. 
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6. Canonical Transformations 

Canonical transformations {x.p)^ {X{x,p), P{x,p)) preserve the phase-space vol- 
ume (area) element (again, take h= 1) through a trivial Jacobian, 

dXdP = dxdp {X, P} , (54) 

ie, they preserve Poisson Brackets 

du dv du dv 
dx dp dp dx 



{X,P},p = l, {x,p}xp = l, {x,p} = {X,P}. (56) 

Upon quantization, the c-number function hamiltonian transforms "classically" , 
H{X,P) = H{x,p), like a scalar. Does the ★-product remain invariant under this 
transformation? 

Yes, for linear canonical transformations'^*', but clearly not for general canon- 
ical transformations. Still, things can be put right, by devising general covariant 
transformation rules for the *-product^®: The WF transforms in comportance with 
Dirac's quantum canonical transformation theory-"^. 

In conventional quantum mechanics, for canonical transformations generated by 
F{x,X), 

the energy cigcnfunctions transform in a generalization of the "representation- 
changing" Fourier transform^, 



,{x) = Ne j dX e'^(^'^) ^e{X) . (58) 



(In this expression, the generating function F may contain h corrections to the clas- 
sical one, in general — but for several simple quantum mechanical systems it manages 
not to^°.) Hence^^, there is a transformation functional for WFs, T{x,p; X, P), such 
that 

f{x,p) = jdXdPT{x,p;X,P) Q)T{X,P) = j dXdP T{x,p;X,P) J^{X,P) , (59) 
where 

T{x,p;X,P) (60) 
= dYdy exp (^-lyp + iPY - iF*{x - |, X - |) + iF{x + |, X + |)^ . 

Moreover, it can be shown that^^, 

H{x, p) ★ T{x, p; X, P) = T{x, p; X, P) ®H{X,P). (61) 
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That is, if satisfies a ©-genvalue equation, then / satisfies a ★-genvalue equation 
with the same eigenvalue, and vice versa. This proves useful in constructing WFs 
for simple systems which can be trivialized classically through canonical transfor- 
mations. 



7. Omitted Miscellany 

Phase-space quantization extends in several interesting directions which are not 
covered in such a short introduction. 

Complete sets of WFs which correspond to off-diagonal elements of the den- 
sity matrix (sec next section), and which thus enable investigation of interference 
phenomena and the formulation of quantum mechanical perturbation theory are 
covered in 47,5i,40 rpj^g spectral theory of WFs is discussed in ii'52,40 gpjj^ jg 
treated in ref . 

Mathematical, uniqueness, and computational features of ★-products in various 
geometries are discussed in ii.54,55,56,57,58,59,60 

Alternate phase-space distributions, but always equivalent to the WF, are re- 
viewed in 9.61,55,32_ 



8. The Weyl Correspondence 

This section summarizes the bridge and equivalence of phase-space quantization to 
the conventional formulation of Quantum Mechanics in Hilbert space. 

WeyP introduced an association rule mapping invertibly c-number phase-space 
functions g{x,p) (called classical kernels) to operators in a given ordering pre- 
scription. Specifically, p i— > p, a; i— > y, and, in general, 

'S(y,p) = / d'T'i'^dxdp g{x,p)exp{iT{p-p)+ia{-^-x)). (62) 

The eponymous ordering prescription requires that an arbitrary operator, regarded 

as a power scries in y and p, be first ordered in a completely symmetrized expression 
in y and p, by use of Heisenberg's commutation relations, [y, p] = ih. A term with 
m powers of p and n powers of y will be obtained from the coefficient of r^cr" in 
the expansion of (rp -|- cry)"'+". It is evident how the map yields a Weyl-ordered 
operator from a polynomial classical kernel. It includes every possible ordering with 
multiplicity one, eg, 

e/x^ ^ p2y2 + y^p^ + pypy + pfp + ypyp -I- yp^y . (63) 

In this correspondence scheme, then, 

Tr6 = j dxdp g. (64) 
Conversely^, the c-number phase-space kernels g{x,p) of Weyl-ordered operators 
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©(y, p) are specified by p i— > p, y i— > a;, or, more precisely, 

g{x,p) = j drda e'(^P+'^^)Tr (e'^^^P+'^f)©) 

= ^Jdye-'y''{x + ^y\(&{T:,p)\x-^y), (65) 
since the above trace reduces to 

j dz e'^<^^/2(2|e-''^'=e-'^P©|0) = '^t^ j dz{z - /ir|©|^)e''"(^'^/2-z)_ (gg) 

Thus, the density matrix inserted in this expression^ yields a hermitean gener- 
alization of the Wigner function: 

fmn{x,p) = ^ j dy e~''yP{x + ^y\il)n){ipm\x - ^y) 

= dye-'^^rmix - ^y)Mx + |y) = /:„(x,p) , (67) 

where the ipmix)s arc (ortho-)normalizcd solutions of a Schrodingcr problem. (Wigner^ 
mainly considered the diagonal elements of the density matrix (pure states), denoted 
above as fm = fmm-) As a consequence, matrix elements of operators, ie, traces of 
them with the density matrix, are produced through mere phase-space integrals^, 

(■^mieiV'n) = y" dxdp g{x,p)fmn{x,p), (68) 

and thus expectation values follow for m = n, as utilized throughout in this review. 
Hence, 

[ipm] expi(c7r + Tp)\ipm) = j dxdp fm{x,p) expi{ax + rp), (69) 

the moment-generating functional^ of the Wigner distribution. 

Products of Weyl-ordered operators are easily reordered into Weyl-ordered oper- 
ators through the degenerate Campbell-Baker-Hausdorff identity. In a study of the 
uniqueness of the Schrodingcr representation, von Neumann''*^ adumbrated the com- 
position rule of classical kernels in such operator products, appreciating that Weyl's 
correspondence was in fact a homomorphism. (Effectively, he arrived at a convolu- 
tion representation of the star product.) Finally, Groenewold'' neatly worked out 
in detail how the classical kernels / and g of two operators 3" and C5 must compose 
to yield the classical kernel of 

= / d^dndi,'dn'dx'dx"dp'dp"f{x',p')g{x",p") 

X expi(e(p - p') + r?(y - x')) expi(^'(p - p") + r?'(y - x")) = 
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{2n) 



J d^dnd£,'dn'dx'dx"dp'dp"f{x',p')g{x",p") exp j ((c + 0\> + {'n + v')t) 



X expi( -ip' - r]x' - Cp" - rj'x" + ^{^t]' - t]C) ) • (70) 



Changing integration variables to 

2 2 2 2 

^' = -{x-x'), i = T - -{x-x'), r]' = -{p'-p), r] = a- -{p' -p), (71) 

reduces the above integral to 

= j drdadxdp expi (t(p - p) + <7(r - x)) {f-kg){x,p), (72) 

where f -kg is the expression (11). □ 
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Appendix A Selected Publications on Phase-space Quantization 

Implicitly, the bulk of the formulation is contained in Groenewold's and Moyal's 
seminal papers. But this has been a slow story of emerging connections and chains 
of ever sharper reformulations, and confidence in the autonomy of the formula- 
tion arose very slowly. As a result, attribution of critical milestones cannot avoid 
subjectivity: it cannot automatically highlight merely the earliest occurrence of a 
construct, unless that has also been compelling enough to yield an "indefinite stay 
against confusion" about the logical structure of the formulation. 

H Wcyl (1927)'^ introduces the correspondence of "Wcyl-ordcrcd" operators to 
phase-space (c-mimber) kernel functions (as well as discrete QM application of 
Sylvester's (1883)'^'^ clock and shift matrices). 

J von Neumann (1931)^^, in a technical aside off a study of the uniqueness of 
Schrodinger's representation, includes an implicit version of the ^-product which 
promotes Weyl's correspondence rule to full isomorphism between Weyl-ordered 
operator multiplication and ★-convolution of kernel functions. 

E Wigner (1932)^ (and L Szilard, impublishcd) introduces the eponymous phase- 
space distribution function controlling quantum mechanical diffusive fiow in phase 
space. It specifies the time evolution of this function and applies it to quantum 
statistical mechanics. 

H Grocncwold (1946)^. A seminal but somewhat imapprcciatcd paper which 
fully understands the Weyl correspondence and produces the WF as the classical 
kernel of the density matrix. It reinvents and streamlines von Neumann's construct 
into the standard ★-product, in a systematic exploration of the isomorphism between 
Wcyl-ordercd operator prodiicts and their kernel function compositions. It further 
works out the harmonic oscillator WF. 

J Moyal (1949)^ amounts to a grand synthesis: It establishes an independent 
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formulation of quantum mechanics in phase space. It systematically studies all 
expectation values of Weyl-ordered operators, and identifies the Fourier transform 
of their moment-generating function (their characteristic function) to the Wigner 
Function. It further interprets the subtlety of the "negative probability" formalism 

and reconciles it with the uncertainty principle. Not least, it recasts the time 
evolution of the Wigner Function through a deformation of the Poisson Bracket into 
the Moyal Bracket (the commutator of ★-products, i.e., the Weyl correspondent of 
the Heisenberg commutator), and thus opens up the way for a systematic study 
of the scmiclassical limit. Before publication, Dirac has already been impressed 
by this work, contrasting it to his own ideas on functional integration, in Bohr's 
Festschrift*^^. 

T Takabayasi (1954)^ investigates the fundamental projective normalization con- 
dition for pure state Wigner functions, and exploits Groenewold's link to the conven- 
tional density matrix formulation. It further illuminates the diffusion of wavepack- 
ets. 

G Baker (1958)^ envisions the logical autonomy of the formulation, based on pos- 
tulating the projective normalization condition. It resolves measurement subtleties 
in the correspondence principle and appreciates the significance of the anticommu- 
tator of the ★-product as well, thus shifting emphasis to the ★-product itself, over 
and above its commutator. 

D Fairlie (1964)^ (also see ^^) explores the time-independent counterpart to 
Moyal's evolution equation, which involves the ★-product, beyond mere Moyal 
Bracket equations, and derives (instead of postulating) the projective orthonor- 
mality conditions for the resulting Wigner functions. These now allow for a unique 
and full solution of the quantum system, in principle (without any reference to the 
conventional Hilbert-space formulation). Thus, autonomy of the formulation is fully 
recognized. 

M Berry (1977)^° elucidates the subtleties of the scmiclassical limit, ergodicity, 
integrability, and the singularity structure of Wigner function evolution. 

F Baycn, M Flato, C Fronsdal, A Lichnerowicz, and D Sternheimer (1978) 
analyzes systematically the deformation structure and the uniqueness of the formu- 
lation, with special emphasis on spectral theory, and consolidates it mathematically. 
It provides explicit solutions to standard problems and introduces influential tech- 
nical innovations, such as the ★-exponential. 

T Curtright, D Fairlie, and C Zachos (1998)'^'^ demonstrates more directly the 
equivalence of the time-independent ★-genvalue problem to the Hilbert space for- 
mulation, and hence its logical autonomy; formulates Darboux isospectral systems 
in phase space; establishes the covariant transformation rule for general nonlinear 
canonical transformations (with reliance on the classic work of P Dirac (1933)^; and 
thus furnishes explicit solutions of nontrivial practical problems on first principles, 
without recourse to the Hilbert space formulation. Efficient techniques, e.g. for 
perturbation theory, are based on generating functions for complete sets of Wigner 
functions in T Curtright, T Uematsu, and C Zachos (2001)^°. A self-contained 
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derivation of the uncertainty principle is provided in T Curtright and C Zachos 
(2001)37. 

M Hug, C Menke, and W Schleich (1998)*'' introduce techniques for numerical 
solution of ★-equations on a basis of Chebyshev polynomials. 
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